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Generation of entanglement between two qubits by scattering an entanglement mediator is dis- 
cussed. The mediator bounces between the two qubits and exhibits a resonant scattering. It is 
clarified how the degree of the entanglement is enhanced by the constructive interference of such 
bouncing processes. Maximally entangled states are available via adjusting the incident momentum 
of the mediator or the distance between the two qubits, but their fine tunings are not necessarily 
required to gain highly entangled states and a robust generation of entanglement is possible. 
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I. INTRODUCTION 

Entanglement plays an essential role in the ideas of 
quantum information, like quantum computation, quan- 
tum communication, quantum cryptography, and so on 
0. Efficient generation of entanglement hence consti- 
tutes an essential element for the realization of such ideas, 
and various schemes have been proposed. Entanglement 
vifould be most simply and naturally generated by a di- 
rect interaction between two entities carrying quantum 
information i.e. between two qubits. There are, how- 
ever, several setups in which the two qubits are separated 
from each other, beyond the range of the direct interac- 
tion, but an entanglement is to be shared between them. 
In such "mediator" would be convenient to en- 

tangle the separated qubits, and several schemes have 
been investigated theoretically [111 HH ItIII Ulioll and 
experimentally [Tllll^. 

Among these proposals, we here concentrate ourselves 
on a scheme based on the "successive" interactions of a 
mediator with two qubits to be entangled [E IE S B S 
0,0,0 . Mediator X prepared in a specific state is sent 
to interact successively one by one with qubits A and B, 
each of which is prepared in an appropriate initial state, 
and then the state of X is measured. If X is found in a 
specific state after the interactions, we end up with an 
entanglement between A and B. 

It has been somewhat standard in such schemes to 
assume that mediator X interacts with qubit A(B) for 
a certain time duration ta(b) [E IE IE IM IE 01 1 and 
the generation of a highly entangled state with a high 
probability is accomplished by tuning these times, given 
the coupling constants. It is, however, a subtle problem 
whether the notion of the "interaction time" makes sense 
or not. It might be valid when the size of a wave packet is 
small enough compared with the interaction region, but 
it might not be otherwise. Even in the former case, no 
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rigorous proof has been given explicitly so far, to the best 
of the present authors' knowledge. 

A rigorous approach to this issue would be to formu- 
late it as a scattering problem and such attempts have re- 
cently appeared |E E3 ■ Qubits A and B are shown to be 
entangled after scattering mediator X. In this approach, 
the notion of the interaction time is not necessary; actu- 
ally, such a time is either absent in the stationary-state 
treatments or automatically given by the physical situa- 
tion under consideration in the wave-packet scatterings 
^El- At this point, it is worth noting that X can be re- 
flected by each of the qubits with a certain probability. 
This effect is not taken into account in the standard for- 
mulations with the interaction times. The wave reflected 
by B is partially directed to A and some of its portion is 
reflected back to B. Such a sequence of bounces between 
A and B gives rise to a resonant scattering of mediator 
X. 

The purpose of the present work is to clarify the ef- 
fects of the ID resonant scattering on the entanglement 
generation 0| : resonance can enhance the degree of en- 
tanglement. It is revealed that the distance d between the 
two qubits is an important parameter to gain an entan- 
glement efRciently, and maximally entangled states are 
available with finite probabilities by adjusting d or the 
incident momentum k of X, in a wide parameter region 
of the coupling constants. Fine tunings of k or d, how- 
ever, are not necessarily required and robust generation 
of highly entangled states is possible. 

II. SCATTERING A MONOCHROMATIC WAVE 

We consider the following Hamiltonian in ID space: 

u , I (X) (A) , (X) {A)s.. , 

Im 

+ gB(Tf Vi^) + o^''of)b{x - d/2). (1) 

Qubits A and B are placed at a; = —d/2 and d/2, re- 
spectively, and the quantum information is encoded in 
their spin states |T)a(b) a-ud |i)A(B)i which are flipped 
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FIG. 1: Mediator X prepared in ||)x is scattered by qubits A 
and B placed at a: = —d/2 and d/2, respectively, and initial- 
ized in |tt)AB, and |t)x is to be detected by either of the two 
detectors on both sides, with a certain probability. 



by the ladder operators a^^^^\ A and B are initially 
prepared in |tT)AB and mediator X polarized in |J,)x is 
injected from the left. X, whose position and momentum 
are represented by the operators x and p, respectively, 
propagates according to the Hamiltonian ^ and is scat- 
tered by A and B (Fig. Since the interactions be- 
tween mediator X and qubits A, B preserve the number 
of spins in the up state |t) among A, B, and X, if X after 
the scattering is detected by either of the two detectors 
on both sides and found in |t)x, either state of A or B 
is flipped down and an entangled state is generated [see 
(|13|l below] . Note that we are looking at an extreme case 
where the potential barriers are much thinner than the 
extension of the wave packet of X. 

Let us find a scattering state of the Hamiltonian 



H\E)=E\E), 



(2) 



which would describe the scattering of a wcU-mono- 
chromatized incident particle X. We introduce the wave 
functions 

= (x(.Ti|®AB(TTl) \E) = (xUniE), etc. (3) 
The solution is given in the form 

(j^ik{x+d/2) ^_ ^^-ik(x+d/2) < -d/2), 

u{x) = <^ (4) 

for each spin state, where k is the momentum (apart from 
h) of X and = tv^k"^ /2m. 

We are interested in the entanglement generation from 
the initialized state |tT)AB by injecting X in the spin 
state |i,)x- In such a case, only the three wave func- 
tions U|;||(a;), M|;|x(a;), and it|;||(a:) arc involved in the 
problem. Imposing the continuity condition on the wave 
functions, the scattering problem Q is solved under the 
boundary conditions 

/X;TT = N, = 0, /T;XT = (5) 

to yield the transmission and reflection amplitudes 
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where 
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and riA(B) = ™5A(B)/^^fc- Note that X does not feel 
A(B) when they are in the same spin states. The coeffi- 
cients iA(B) and ?'a(b) respectively describe the transmis- 
sion through and the reflection from A(B) without spin 
flip and /a(b) the transmission/reflection with spin flips, 
when X is in a different spin state from A(B). 

It is interesting to expand the amplitudes in ll^Jl as 



power series m e 
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etc., where e represents the phase factor gained for the 
rightward/leftward propagation of X between A and B, 
and each term of the expansions reveals how X goes back 
and forth between A and B, before flipping its spin by 
the interaction with A or B. 

It is also worth looking at the coefficients iA(B)5 ?'A(B)j 
and /a(b) for the single potential A(B) and comparing 
the present approach with the ordinary formulation with 
the interaction times 0, IE S S IS • When we con- 
centrate ourselves on a transmitted particle (like in the 
approach with the interaction time), a state of X and 
A, e.g. II; t)xA, is "rotated" after the transmission of X 
through potential A like 

U;T>xA^<AU;T>xA + ./A|T;i)xA, (9) 

apart from the normalization. The counterpart in the 
approach with the interaction time reads ^ |8j 

U; T)xA cos^ataU; T)xa - isingATAlT; i)xA, (10) 

and there exists the connection between the two formu- 
lations: 

1 . ^A 



C0S5ATA 



sm^A^A 



(11) 
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These relations show that qata ranges from to 
7r/2 and a higher/lower momentum corresponds to a 
shorter/longer interaction time. 

In practice, however, one should also take account of 
the probability for such event, transmission, to occur. To 
achive the complete flip smgATA = 1, for instance, the 
correspondence (|llf) suggests that 17a — > oo is required. 
But the probability of the transmission 



\tA? + \lA? 



1 



(12) 



vanishes in this limit, meaning that the complete flip is 
not possible. A larger rotation angle ^a^a requires a 
larger Oa but it is available with a smaller probability. 
This restricts the applicability of the formulation with 
the interaction times, at least for the delta-shaped po- 
tential. 



III. GENERATION OF ENTANGLEMENT 

Now we are ready to discuss the entanglement gener- 
ation. If a transmitted/reflected particle in |t)x is de- 
tected by the detector on the right/left, an entangled 
state 



|vI^t)AB«rT;Tilti)AB+TTaT|iT) 
(for the former case) or 

|*,)ABCxi?T;TilTi)AB+i?T;iTUt) 



AB 



AB 



(13a) 
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(for the latter) is established. In the present model, the 
concurrences Il5l of these states are the same, 
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(14a) 



l + 4fi|(l + r2^)sin^fcd, 
(14b) 



and so are the probabilities of the generations of the en- 
tangled states (|13a|l and Ijl3b|) . namely the probability 
of detecting a transmitted particle in the |t)x state and 
that of detecting a reflected one. 



P 



I^T;TiP 
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2 I An2 o2 l^ I n2 ^ „;„2; 



% + 4r2^f7^(l + %) sin'^kd 



(15) 

Note that the probability of either of the two detectors 
detecting X in the state |t)x is given by 2P, and therefore 
the maximal value of P is 1/2. 

Notice here that the distance between the two qubits, 
d, enters the formulas through the exponential factor 
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FIG. 2: (a) The concurrence C given in 11411 and (b) the 
probability P given in 11511 as functions of k and qa when 

qb ~ 3 [h^Tv/md]. 



^2ikd^ which is responsible for the resonant scattering. 
See Fig. |21 where the resonance structures are observed 
in the momentum dependences of the concurrence C and 
the probability P. By adjusting the incident momentum 
k, one can generate a highly entangled state with a finite 
probability. At the same time, fine tunings of param- 
eters are not necessarily required. For example, let us 
look at a case where ^a — SB- Figure Ola) shows that 
C and P do not oscillate strongly for large k and a high 
concurrence is available in a wide range of the incident 
momentum k. This is because the oscillating factor e^'*^'' 
represents the bouncing of X between A and B and it 
always accompanies the reflection coefficients ta and re , 
which are reduced for large momentum k and suppress 
the oscillations. 

Let us look more closely at the formulas and see how 
a highly entangled state is acquired by the resonant tun- 
neling. Figure Olb) shows that the scries expansions in 
© converge very quickly and the first two contributions, 
i.e. the no- and one-bounce processes, are important. For 
9 A = 9b, the no-bounce process, i.e. the first terms in ||SJ), 
yields |T|;|J < \T^,i^\ (|i?T;Til < l-RT;iTl)- If momen- 
tum k is adjusted so as to satisfy kd = {u = 1,2, . . .), 
the one-bounce process, i.e. the second terms in in- 
terferes with the no-bounce process constructively for the 
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FIG. 3; (a) The concurrence C given in H14^ (solid) and the 
probabihty P given in 1151 (dashed) as functions of the inci- 
dent momentum k when (?a = 3b = 3 [h^Tv/md]. (b) The same 
as (a) but only the contributions up to the (n + l)th term in 
eacli of the series expansions (|HJ| are taken into account, where 
n is the number of bounces between A and B. 



It; ID component but destructively for the |t;iT) com- 
ponent (since iA(B) > while ?'a(b) < 0). As a result, the 
difference between |r|;-f J and \T-^;i^\ (|i?T;Til and |i?T;iTl) 
is reduced by the interference and the concurrence is 
made close to unity. In fact, the ratio a introduced in 
()14|l is evaluated up to the one-bounce contribution to 
hea = l + Sa=l + + 2^"^ + 2^1'^) under the 

conditions gx = 3b and kd ^ vn (v ~ 1,2, . . .), and the 
deviation of the concurrence from unity is about ((5a)^/2, 
which is very small for Q (= = ^^b) ^ 1 (almost per- 
fect compensation). Other processes with more bounces 
complete the exact unit concurrence. Note that, in Fig. 
13 where = 5b = 3 [^i^tt/toc?], the condition Q < 1 
corresponds to fc > 3 [-/r/d]. 



We first look at the probability H15() . It is possible to 
show that it is optimal for a given pair (fiAji^B) when 
the resonance condition 



sin^fcd 



(16) 



is satisfied. The probability P under this condition is 
plotted in Fig.QJa) as a function of flA and rie, together 
with the concurrence C. The maximum value of the prob- 
ability P is given by 



at (r2A,f^B) = (l/%/2,oo), (17) 



while the unit concurrence C = 1 under the resonance 
condition (|16|l is achieved when 



1 



2ni' 



(18) 



If one does not stick to a high probability, the concur- 
rence C is optimized as follows. The concurrence (|14f) 
takes its maximum C = 1 at a = 1 and decreases mono- 
tonically as a leaves this optimal point. Since 



2i<«<2^(i + 2si|, 



(19) 



for a given pair (r^Aii^B), the unit concurrence C 
(i.e. a = 1) requires 



^ < 1 < ^(l + 2f}|), sin^fcd: 
Ub "b 



"b "a 



(20) 

See Fig. E^b) , where the region for the unit concurrence 
C = 1 is shown, together with the corresponding prob- 
ability P. In particular, the unit concurrence C = 1 is 
always available when = 5b, by adjusting the incident 
momentum k or the distance between the two qubits, d. 
The optimal probability P for the unit concurrence C = 1 
is realized at the point indicated by a dot in Fig. ^b) , 
which is evaluated to be 
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IV. OPTIMIZATION 



ANOTHER MODEL 



Let us next survey the whole parameter space and 
discuss how to optimize the concurrence C in H14(l 
and the probability P in (|15|) . There are essentially 
three independent (dimensionless) parameters r2A(B) = 
m5A(B)/^^fc and kd. 



Similar analyses are possible for the Hamiltonian 

2 

H^^+ 5Acr(^) • (T<-^^6ix + d/2) 
2m 



+ 5BO-(^) • cr(^)j(a; - (i/2). 
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FIG. 4: (a) Contour plots of the probability P (dashed) given in 11511 and the concurrence C (solid) given in 1141 as functions 
of S7a and S^b under the resonance condition sin^fcd = 1. The probability P takes its maximum value Pmax = 1/2 at 
(r^A, Ob) = (l/\/2, oo), and the concurrence is unity C = 1 on the line SIa ~ f2B/(l + 2S7|). (b) Contour plots of the optimal 
concurrence C for given (SIajS^b) (solid) and the corresponding probability P (dashed) as functions of Oa and Qb- The unit 
concurrence C = 1 is available in the gray region f2B/(l + 2r2|) < f^A < f^B- In the left region Oa < ^•ajiX + 2r2|), the 
concurrence C is optimal when sin^fcd = 1, and in the right region JIa > f^B, it is optimal when sin^fcd — 0. The optimal 
probability P for the unit concurrence C = 1 is realized at the point indicated by a dot, {Q.p^,Q.^') ~ (0.33, 1.07), and is given 
by Popt ==; 0.37. 



This type of interaction is considered in Ref. |^ 0| . 
this case, the amplitudes are given by 
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In contrast to the previous model, X feels potential A(B) 
even when X is in the same spin state as A(B). t'x[B) ^^^^ 
^A(B) coefficients for the transmission through 

and the reflection from the single potential A(B) with- 
out spin flip when X is in the same spin state as A(B), 
iA(B) and ''A(B) are those when X is in a different spin 
state from A(B), and /a(b) describes the transmission 
and the reflection with spin flip when X is in a different 
spin state from A(B). For this model, the concurrence Ct 
and the probability Pt by detecting a transmitted par- 
ticle in ||)x on the right are different from Ci- and Pi- 
by detecting a reflected one on the left. Again, the de- 
nominators in H23|l and H25|l represent bouncing of X be- 
tween A and B (which becomes clear by expanding them 
as power series in e^^^'^) and give rise to the resonant 
tunneling. Similar resonance structures to the previous 
example are observed in the concurrences and the prob- 
abilities, as shown in Figs. O and |H| The oscillations of 
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FIG. 5: The concurrences (above: a,c) and the probabilities (below; b,d) for the Hamiltonian 12211 as functions of k and qa 
with qb = 1.5 [h^Tv/md], when a reflected particle is detected (left: a,b) and when a transmitted particle is detected (right; 
c,d). 
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FIG. 6; The concurrences (solid) and the probabilities (dashed) for the Hamiltonian H22|l as functions of the incident momentum 
k with flA = 3b = 1.5 [fi,^7r/md], when a reflected particle is detected (a) and when a transmitted particle is detected (b). 



the concurrences and the probabilities arc reduced for a 
large momentum k and a robust entanglement generation 
is available (Fig. ^ , due to the same mechanism as in the 
previous example. 



VI. SUMMARY 



In this article, we have investigated the entanglement 
generation by the resonant scattering. The resonance 
effects are clarified and the optimization of the entan- 
glement generation is discussed. The interference of the 
bouncing and non-bouncing processes can enhance entan- 
glement. Maximally entangled states are available with 
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finite probabilities in a wide parameter region. The de- 
gree of entanglement is optimized by adjusting the mo- 
mentum of the mediator, but its fine tuning is not nec- 
essarily required. 

One of the interesting future subjects would be to clar- 
ify the effect of the size of a wave packet on the generation 
of entanglement. Scattering of plane waves is discussed 
in this article, but the models are solvable also for wave 
packets 0]. It is possible to discuss the resonant scat- 
tering fully dynamically. Introduction of the width of the 
potential would also be interesting to explore the validity 
of the ordinary formulation with the interaction times. 
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